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SUMMARY 


During the re-entry phase of a manned satellite, some equipment for 
continuous on-hoard indication of position will he required. Since radio 
and radar may he useless during part of the re-entry, inertial guidance 
equipment may he required. Such equipment, however, has an inherent 
instability in the computation of altitude . 

The present study of an inertial guidance system shows that for 
reasonable values of initial- condition errors and equipment biases, the 
res ul tant position Indication errors will not become excessive unless 
the re-entry maneuver time is greater than 45 minutes to an hour. Further, 
the position indication error caused by accelerometer uncertainties can 
be reduced by switching off the accelerometers until their output becomes 
significantly greater than their uncertainty. 


INTRODUCTION 


For a mann ed satellite vehicle re-entering the earth's atmosphere 
to lan d at a particular spot on the earth, the pilot must continuously 
know his spatial position, velocity, and attitude. Before re-entry, 
this information can be obtained either from on-board sources or, in the 
case of position and velocity, from ground-based radars with telemetering 
to the vehicle . 

However, during the actual re-entry, the vehicle will be surrounded 
by a layer of hot ionized gas and, at altitudes less than about 350,000 
feet, may not be able to communicate with the ground via radio or radar. 
Radio contact with the ground probably would be re-established by the 
time the vehicle reached an altitude of about 150,000 feet. Thus, for 
the interval between 350,000 and 100,000 feet, a self-contained on-board 
source of position information - an inertial navigation system - probably 
will be required. A system consisting of a gyro reference platform to 
maintain an orientation in space, an accelerometer package mounted with 
the gyros to measure the vehicle accelerations, and a computer to deter- 
mine position, velocity, and attitude from the outputs of the accelero- 
meters and gyros would be capable of measuring position and velocity 



fairly accurately over long periods of time if it had a suitable means 
of determining altitude. Present radar or barometer systems used in 
aircraft for sensing altitude would not be suitable . 

Under the above conditions , the inertial guidance system of a 
re-entry vehicle would have to compute its own altitude. The equations 
for doing so have, at suborbital speeds, an unstable solution, that is, 
one in which any incorrect initial condition or error in acceleration 
measurement will cause an error in altitude indication which will grow 
exponentially with time (see appendix) . The time constant of the solu- 
tion is considerably less than the satellite re-entry time. For this 
reason it seems imperative to make an investigation of the errors that 
may occur in an inertial guidance system during re-entry. In particular, 
can the error caused by incorrect measurement of the acceleration be 
reduced, perhaps by switching out the accelerometers for a period of time 
during the re-entry? The present report presents the results of such an 
investigation . 


LIST OF SYMBOLS 


A 

C D 

CpA 

m 

D 

F 

f 

g 

Sr 


h 

L 

m 

R 

R r 


reference area for aerodynamic force coefficients 
drag coefficient 

vehicle constant relating drag to airspeed and air 
density 

aerodynamic drag per unit mass, ft/sec 2 
force per unit mass, ft/sec 2 

nongravitational force per unit mass , ft/ sec 2 

gravitational attraction of earth per unit mass, ft/sec 2 

gravitational attraction of earth per unit mass at sea 
level, 32.1279 ft/sec 2 

altitude above sea level, ft 

aerodynamic lift per unit mass, ft/sec 2 

mass of vehicle, slugs 

distance from satellite to center of earth, ft 
radius of earth, 0.209102x10 ft 
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a, 3,7 
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oa 
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OU 


thrust per unit mass from retrorocket, ft/sec 2 
time, sec 

uncertainty of accelerometer reading, ft/sec 2 

velocity of satellite relative to the air mass, ft/sec 

axes of satellite coordinate system 

axes of indicated satellite coordinate system 

axes of platform coordinate system 

axes of platform coordinate system at t = 0 

unit vectors along satellite coordinates 

velocity increment at initial altitude 

flight-path angle - angle from true horizontal to 
V^, positive up, radians 

air density, slugs/cu ft 

supplement of angle between Xp and T, radians 

range angle, measured at the center of the earth 
between a reference line and the vertical, radians 

specific times during a re-entry^ sec 

angular velocity with respect to inertial space, 
radians/ sec 

gyro drift rate, radians/sec 

earth angular velocity, 0 . 7292 x 10 4 radian/sec 
Superscripts 

vector 

time derivative 



Subscripts 


i indicated quantity 

0 value at t = 0 

e error between indicated and true 

a components of vector along a and y axes 
x,z components of vector along X and Z axes 

m measured quantity 

GENERAL CONSIDERATIONS 


The approach used was to solve both the equations of actual motion 
of the satellite in response to the actual forces and the equations of 
indicated motion (used by the inertial guidance equipment) in response 
to measured forces . The latter equations have in them the initial- 
condition errors and equipment biases that might reasonably be expected 
to occur. The two sets of results were then compared to ascertain the 
position and velocity errors of the inertial guidance equipment. 

These errors were studied as functions of initial- condition errors, 
equipment biases, orbiting altitude, vehicle lift-to-drag ratio, and 
re-entry maneuver. The effect of having the accelerometers switched out 
for various periods during re-entry was also ascertained. To simplify 
the equations , it was assumed that the vehicle moves easterly in the 
equatorial plane and that the earth is spherical. These simplifying 
assumptions are justified since we are trying to isolate the errors 
caused by the inertial navigation system. Similarly, small discrepancies 
between the actual and assumed atmospheric density profiles will not 
affect the errors caused by the inertial navigation equipment . 

A set of basic conditions was chosen for use during this study and 
then various parameters were changed to study their effects on the errors . 
The standard vehicle was assumed to have an L/D of 0.5 with a C^A/m 

of O.533. The "basic initial altitude was 677,000 feet and the basic 
velocity increment, AV, selected was 250 feet per second. This and all 
other values of AV were developed by applying a lg thrust at an angle 
of L5 0 down and to the rear of the flight path, for that length of time 
which was necessary to develop the desired AV . 

For these conditions, the position errors resulting from accelerometer 
biases and gyro drift, both singly and collectively, were determined by 
computation, on a digital computer. The effects of initial-condition 
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errors were also determined, both singly and collectively, and finally 
the effect of switching out the accelerometers for various periods of 
time immediately following retrothrust was determined. 

The effects of accelerometer biases and gyro drifts were then 
determined for two other altitudes, two other vehicles (a zero-lift 
vehicle with CpA/m of 1.0, and a glider with CpA/m of O.I 58 and 
L/D of l.o), and two other values of AV (at the same angle). For all 
of these digital-computer runs the effects of switching out the 
accelerometers were determined . 

As used in this report, uncertainty of a quantity implies an 
imperfect knowledge of that quantity. Thus, for example, the acceler- 
ometer output is equal to the actual value of the acceleration plus the 
accelerometer uncertainty. An uncertainty has a random component plus 
a quasi-steady component. The latter component is referred to in this 
report as bias . 


Coordinate Systems 


Two basic coordinate systems are used in this report. Both are 
right-handed, orthogonal, and centered in the satellite, with one axis 
forward and one axis down. One is the normal aircraft system (here 
called the satellite coordinate system) with coordinates a, 3, and 7 , 
in which the 7 axis is along the true vertical (toward the center of 
the assumed spherical earth) and the a axis is forward. The other 
coordinate system, called the indicated satellite coordinate system, is 
derived by the on-board computer, and is used to determine the indicated 
position of the vehicle. It has coordinates X, Y, and Z, with the Z 
axis along the indicated vertical. Since the problem is solved in a 
single plane, the Y and (3 axes coincide. 

Two additional coordinate systems are used as references . One is 
the platform coordinate system, with coordinates Xp, Yp, and Zp, such 
that Yp ' coincides with Y and 3 . The platform is designed to be 
inertially fixed, and the on-board computer assumes that it is . However, 
there will actually be some rotation of the platform caused by gyro drift. 
The other reference coordinate system represents the location of the 
platform coordinate system at time zero and is inertially fixed. Its 
coordinates are X po , Y po , and Zp 0 * Again, Y po coincides with Y, (3, 
and Y p . The angular motion of the platfoim system with respect to the 
inertially fixed system is about the Y p axis and represents gyro drift. 
The Zpo axis is chosen along the initial true vertical. 

Figure 1 shows the four coordinate systems and the various vectors 
used in this report. The vectors (3, Y, Y p , Y po , and w are perpendicular 
to, and coming out of, the page. 
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DERIVATION OF EQUATIONS 
Equations of Motion of the Satellite 


The motion of the center of mass of the satellite can he expressed 
in an inert ially fixed coordinate system as: 

f = F (1) 

where F is the total force per unit mass exerted on the body. This 
force can he divided into two parts: g, the gravitational attraction of 

the earth, and f which consists of the lift, L, drag, D, and the thrust, A 

T. All forces other than those mentioned above have been neglected. 

Equation (l) can he written in an earth-centered coordinate system, 
rotating at an angular velocity with respect to inertial space of u by 
means of the rule of Coriolis, as follows (ref. 1, p. 3^2) : 


R + wx(coxR)+2cjxR + u)xR=F (2) 

where 


F = g + f = g + L + D+ T 


This equation is also valid when expressed in terms of a satellite 
coordinate system, where y is toward the earth center, and w Is the 
orbital angular velocity. 

Expressing the quantities used in equation (2) in this satellite 
coordinate system gives: 


* * 

R = -(R r + h) 7 R = -h 7 


R = -h/ 


CO = “03 



g = 


(i“Thp 


7 


( 3 ) 


It should be noted that the angle 0 is defined so that it will 
increase positively as the satellite moves around the earth, and does 
not follow the usual right-hand rule for angles and angular velocities 
in the satellite coordinate system. Also, 0 is measured with respect 
to an inertially fixed reference, and not with respect to the rotating 
earth . 
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Substituting equation (3) into equation (2), we get 


-I17 + (R r + h) 9 7 + 2hda + (R r + h)0ct = 


Sr^r 

(H r + h Y 


7 + f a a + f 7 7 


( 4 ) 


where f and t y are the components of f along a and 7, respectively , 
and fp = 0, since we are assuming a two-dimensional problem. 


These equations can now be separated into their a and 7 components, 
giving 


h - (R r + h) 9 - - T-— C . 2 " f 7 

(R r + h) 

(R r + h)e + 2h0 = f a 


( 5 ) 


which describe the actual motion of the center of mass of the satellite 
in response to the forces g, f a , and f^^ 


The force f can be determined by summing vectorially the lift, 
drag, and thrust of the vehicle as shown in equation (2) . 


The equations for the magnitude of D and L are 

D = (l/2)pV A 2 (C D A/m) 

L = (L/D)D 


( 6 ) 


These quantities are directed antiparallel and perpendicular , 
respectively, to the vector representing the velocity of the satellite 
with respect to the air mass. The quantity CjjA/m is a constant of the 

vehicle, and it is assumed that the vehicle angle of attack is controlled 
so that the quantity L/D is also constant. The air density, p, was 
taken from the standard ARDC atmosphere (ref. 2) for altitudes below 
332,000 feet, and from the Sputnik atmosphere for higher altitudes 
(ref. 3). 

The retrothrust, T, is used only during the initial portions of the 
re-entry to change the vehicle motion from orbital to re-entry. It is 
applied along a body axis, which is assumed to be maintained, during the 
thrust period, in a specified orientation, q), with respect to the inertial 
platform. Following the thrust period, the body is controlled to maintain 
constant L/D. 



8 


Resolving f into f a and f ^ gives: 

f_ = -L sin 7 - D cos y - T cos(e + cp - 

O- A A 

f = -L cos 7„ + D sin y + T sin(0 + 9 - w,t) 
7 A A u - 


( 7 ) 


Then substituting equation (6) into equation (7), we get 

f a - - | PV/ ^ (cos 7 a + p sin 7a ) - T cos( 0 + cp - ^t) 

f 7 = \ P V A 2 ( Sln 7 A - i COS 7 a) + T Sln(0 + 9 “ W <i t) 


The following sketch shows ^the derivation of V A and 7 A from 
satellite and earth velocities. 



A 


From the sketch, we get the following equations: 


V a 2 = h 2 + (R r + h ) 2 (0 - “ef 


^ -1 h 

7. = tan — — 7 

^ (R r + h)( 0 - w e ) 


(9) 


^The constant (R r + h)u e is the velocity of the air mass with 
resnect to inertial space for an atmosphere which is assumed to be 
22 dly attached to the earth and to move with it. Since this velocity 
lies in the equatorial plane, its use assumes that V A and V are also 
in this plane and that the vehicle is thus in an equatorial orbit mov ng 

in an easterly direction. 


vo vo oj 


m o\ CT\ 


2U 


9 


Equations (5), (8), and (9), and the atmospheric table describe completely 
(within the limits of the assumptions) the motion of the satellite. 


Inertial Navigation Equations 


A 


Since equations (5) describe the motion of the satellite in response 
to gravity and the force f, these equations can also be used by the on- 
board computer to determine satellite position in response to the measured 
forces and computed gravity. The equations then become, expressed in 
computer quantities, 


^i “ ( R r + ^i)®i ~ zm 

®i( R r + h i) + 2 ^i®i = f xm 


SrRr 2 

(R r + hi ) 2 > 


( 10 ) 


Here f* and f are components of the measured forces, in the 

xm zm r ^ y 

indicated horizontal and vertical directions, respectively* The 
accelerations are actually measured in platform coordinates. The accel- 
erometer outputs can thus he obtained by resolving the true forces, 
f a and f y 9 into platform X and Z components through the angle (-0 + uj^t) 

(see fig. 1) and then adding in the uncertainties of measurement 

U-xp and U Z p. These accelerometer outputs must then be resolved through 

the angle 0. to obtain f and f . Both of these resolutions were 
1 xm zm 

done simultaneously, giving the following equations: 


f 


xm 


f 


zm 


f a cos (-0 + co^t + Qi) + fy sin (-0 + + 9$) 

+ Uxp cos 9± + U 2 p sin 

f^, cos (-0 + + 0^) - f^ sin (-0 + w^t + 0^) 

+ u zp COS e ± - Uxp sin e ± 


k (11) 


J 


2 

The quantities f^ and f zm are the components of the accelerometer 
output (force per unit mass) in indicated horizontal and vertical direc- 
tions; g r R r 2 /(R r + t^) 2 is the computed gravity term which must be added 
to the accelerometer output to determine the actual linear acceleration. 
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Error Equations 


We no-w have equations describing the actual motion of the satellite 
in response to the actual forces (eq. (5))> and the motion as indicated 
by the inertial navigation equipment, in response to the measured forces 
(eq. (10)). The errors of the indicated quantities can be defined as 





J 


(12) A 


Initial Conditions 


Equations which describe the motions of any orbital vehicle (within 
the lim its of our assumptions) have been derived, and now the initial 
conditions must be defined. For this problem, we assume that the vehicle 
lias been traveling in a circular orbit at the specified initial altitude 
h Q (and therefore h Q = h Q = 0) and that the aerodynamic forces f a and f y 
while in orbit are also zero. Using the first of equations (5), we get 



Sr^r 


1/2 


L(R r + h Q rj 


( 13 ) 


The value of 0 O was arbitrarily chosen to be zero. 

The initial conditions for the indicated equations (10) are ideally 
identical to those used for the actual motion of the vehicle. However, 
in practice there will be errors in the indication of position and veloc- 
ity just prior to retrothrust which will alter these initial conditions 
as follows : 


^io - ^o + h eo 


^io _ ^o + k £ o 


®io ®eo 


e lo ■ B o + % 


( 14 ) 
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The effect on the initial conditions of drift of the gyro reference 
coordinate system prior to retrothrust must he considered also, as 
follows: Q is measured from the inertial space reference, Zp 0 , to the 

local vertical, 7; G± is measured from the platform reference, Zp, to the 
indicated vertical, Z. In the presence of gyro drift prior to retrothrust 
Zp and Zp 0 will not he parallel at time zero, and thus, although 
0^ and Q may he equal, the indicated and true verticals will not he 
parallel. Since, in the solution of the problem, Zp Q is assumed along 
Zp at time zero, we can allow for these nonparallel initial verticals 
hy letting he nonzero and equal to the gyro drift angle to he 

considered. 

Thus, G e0 can he used to represent two errors: One is the initial 

error in position indication from either the airborne equipment or ground 
measurements of satellite position. The other is the error in reference 
orientation due to platform drift prior to retrothrust. Both of these 
sources of error are considered. 


METHOD OF SOLUTION 
Computation Equipment Used 


Equations (5), (8), (9), (10), (ll), and (12) must now he solved 
simultaneously to determine the various errors where equations (13) and 
(14) specify the initial conditions. This cannot he done analytically. 
Although the equations could he set up on an analog simulator, the 
answers would he highly inaccurate, since in several places It is neces- 
sary to subtract nearly equal numbers. The use of a digital computer, 
in particular the IBM 704, was indicated as the means of solving these 
equations. Accordingly, these simultaneous ordinary differential equa- 
tions were programmed for the IBM JOk, using single precision floating 
decimal point arithmetic. In the computations, a subroutine applying 
the Adams -Moulton and Runge-Kutta integration techniques was used. 

Double precision arithmetic is found within this subroutine to control 
round-off error. The variable step- size mode of integration was used; 
therefore, the program can determine each interval size sufficiently 
small that the truncation error per integration step is always less than 
some specified value. While the program guarantees that the local error 
Is less than this value, the cumulative error usually exceeds It. For 
this problem, it was desired to use all of the accuracy available from 
the IBM 704. Thus, the truncation error per step was specified as less 
than one part in 10 , so that all the digits of each number were retained. 
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Form of Solution 


The inertial navigation equipment aboard would be used during most 
of the phases of operation of the satellite vehicle. During lift-off 
and boost, the equipment would continuously compute vehicle position, 
using equation (10) . At cutoff, the on-board computer would have deter- 
mined the vehicle velocity and position, and would use these values to 
continuously compute its position in orbit. This could be done using 
equation (10), with f m = f a = 0. However, a more accurate method 
would be to determine the six elements of the orbit and their variation 
with time, and from these to compute the vehicle position. While the 
vehicle is orbiting, this position information would probably be cor- A 

rected by ground measurements of vehicle position, and gyro drift effects 
might be eliminated by stellar monitoring. 

For this study, we have assumed that the on-board computer will use 
equation (10) to determine position and velocity during re-entry. At 
time t = 0, we assume that the computer is solving this equation and, 
since the retrothrust is to begin momentarily, the accelerometer outputs 
f^ and f zm are being used in this equation. 

At time t , arbitrarily chosen to be at 10 seconds, the retrothrust 
rockets are fired. The rocket acceleration is preselected as T, equal 
to 1 g, and a specified value of nominal change in velocity, AV , and 
direction of thrust, <p, is chosen. Then t 2 can be computed as follows: 



The rocket is then turned off at time t 2 . 

For those cases where it was desired to switch out the accelerometer 
input to the computer subsequent to retrothrust firing, a time T g , 10 
seconds later than t 2 , was selected as the accelerometer switch-out 
'time; that is, subsequent to t , the acceleration input to the computer 
was assumed to be zero. The final switch- in of the accelerometers then 
occurred at t 4 . In some cases, this was a preselected time following 
t . In other cases, the outputs of the accelerometers were used as an 
indication of when they should be switched in. That is, the value of 
f and f zm , after resolution into indicated coordinates, was compared 
with the assumed uncertainty of measurement. When f M or f zm became 
significantly greater than this uncertainty, the respective accelerom- 
eter output was applied to the indicated position equations (10) . 

Since f^ and f zm might become significant at different times, the 
X and Z accelerometers were switched in at t 4x and t 4z , respectively. 
These times were not necessarily simultaneous, and either one could 
occur first. 
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For all cases , the indicated position errors are computed down to 
an altitude of 100,000 feet. It is assumed that at this altitude the 
inertial navigation equipment can obtain altitude information from 
barometric pressure and, since communication with the ground has been 
re-established, that position, altitude, and velocity can be checked 
against a radar altimeter and the telemetered ground information, so 
that further indication errors could be corrected. 


Error Magnitudes 


Before computing the indication errors resulting from erroneous 
initial conditions and equipment biases and uncertainties, it is neces- 
sary to know the probable magnitudes of these initial condition and 
equipment errors. 

First, consider initial condition errors. Assume that radar or the 
on-board inertial guidance equipment will compute the burnout position 
to within approximately 200 feet and the corresponding velocity to within 
1 foot per second in magnitude and within 0.00 3° in orientation. Con- 
verting these quantities into the form in which they are used here, we 
have at burnout an altitude error of ±200 feet, a range error of ±10’ 5 
radian, an altitude rate error of ±1 foot per second, and a range rate 
error of ±5x10’ 8 radian per second. 

If the equations previously derived are used to determine the 
vehicle position while orbiting, the indication errors will increase as 
a function of the number of orbits. In particular, after two orbits, 
the range error Q e will increase to ±14x10’ 4 radian and the altitude 
error to around 400 feet while h G and Q e will remain about the same 
size. This assumes that the average effect of drag on the satellite is 
considered in computing the indicated position. If this drag is neg- 
lected, these errors will build up very much larger. However, with 
orbital corrections from the ground, plus a position fix just prior to 
firing of the retrorocket, it is assumed that the values of the initial 
condition errors at t = 0 might be as follows: 

-4 

h €0 = 500 feet 0 €O = 10 radian 

h £0 = 1 foot per second 0 eO = 5x10’ 8 radian per second 

These values were chosen as nominal values for the initial condition 
errors although other values were also studied. 

The effect of equipment uncertainties was also considered. In 
general, the imperfections of gyros show up as gyro drift rate. Rela- 
tively poor gyros were assumed for this study, and a value for of 

2x10“ y radian per second, or about 0.04° per hour, was chosen as a 
nominal drift rate although, again, other values were tried. Imperfec- 
tions in an accelerometer show up as an uncertainty in the output. This 



uncertainty consists of two components, a quasi-steady error, and a 
random error . During any one entry , the system errors due to the random 
component will more or less cancel and can be ignored. The quasi-steady 
error may vary slowly throughout an entry, and will be random over an 
ensemble of entries. Therefore, this term will create a system error 
which must be considered. In order to avoid the necessity for statisti- 
cal measurements of the resultant errors, a specific bias value, constant 
throughout a run, was assumed which would give the largest errors, that 
is, the worst possible case. The results can then be interpreted as 

an envelope of error within which it can be expected that all 
errors will be found. The bias value principally used was 10" g, which 
corresponds to a relatively poor accelerometer. 

It should be noted that because the accelerometer axes are nominally 
fixed in spatial direction, the direction of the bias is also fixed in 
space. However, since the indicated coordinate system is rotating around 
the earth, the biases must be resolved through the indicated range angle 
before being used in the on-board computer, as shown in equation (11) . 

As long as the range angle remains small, the concept of envelope of 
error from accelerometer bias is applicable. However, when the range 
angle becomes greater than 90°, some of the effects of accelerometer 
bias change sign, so that this envelope concept is no longer applicable. 


RESULTS AID DISCUSSION 


Figure 2 shows the flight path profile for each of the various 
conditions. It can be seen in figure 2(a) that different vehicles, 
starting from the same altitude with the same value of AV , follow the 
same path to an altitude of about 300*000 feet before diverging. All 
the trajectories are smooth curves, essentially Keplerian, until this 
altitude is reached. Thus, 300,000 feet is the altitude at which the 
air density becomes great enough to noticeably alter the flight path. 

At altitudes higher than this the flight path is essentially Keplerian, 
with, however, some very small deviations due to the residual air density. 

All these runs were computed down to an altitude of 100,000 feet 
except the Keplerian orbits, and the end of each run is indicated by an 
asterisk. On subsequent figures, an asterisk is placed on each curve at 
a time which corresponds to the vehicle reaching 100,000 feet. 

Figure 2(c) shows the range profile for these same conditions. For 
all vehicles and all entry conditions the Keplerian range profiles are 
essentially identical and are drawn as one curve. The actual range 
profile deviates slightly from this Keplerian value at altitudes below 
300,000 feet. 

Figure 2(d) is a curve of acceleration versus altitude. It shows 
that above 300,000 feet the acceleration level for a given vehicle is 
independent of both the initial altitude and the initial re-entry maneuver. 
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For example, for the standard vehicle, the f a curve is independent of 
AV and h Q , as is the f y curve. Different vehicles, however, have dif- 
ferent acceleration levels. This would indicate that, for a given vehi- 
cle, it would be satisfactory to switch the accelerometers in as a 
function of Keplerian indicated altitude, rather than as a function of 
accelerometer output . 

Re-entry runs were made under many different conditions , with 
different error producing sources of different magnitudes, singly and 
collectively. Study of the resultant errors of these systems showed that 
altitude error seemed to be the best criterion for judging whether or not 
a particular error-producing source caused excessive errors. The actual 
magnitude of the altitude error varied from 0 to ±20,000 feet. 

The effect on the altitude error of various initial- condition errors 
applied individually is shown in figure 3(a) for the standard vehicle 
(CpA/m = O.533, L/D - 0*5) starting from h Q = 677*000 feet with 
AV = 250 feet per second. It can be seen that initial altitude errors 
and range rate errors cause the largest altitude error. Plotted on the 
same figure is a curve for an initial range error of 2.16x10 3 radian. 

This would correspond to a gyro drift rate of 2x10" 7 radian per second 
(0.04°/hr) for a period of approximately 3 hours, or two orbits. The 
resultant altitude error at 100,000 feet is large, indicating that gyro 
drift rates of this magnitude probably can not be tolerated while orbit- 
ing. They must either be corrected by stellar monitoring or by the use 
of better gyros . 

The effects of various other magnitudes of initial-condition errors 
were studied both singly and collectively. In spite of the nonlinear 
equations, the results indicated that the resultant altitude error was 
linearly proportional to the magnitude of the initial-condition error, 
and that the superposition theorem could be used to determine the effects 
of initial- condition errors collectively; that is, the result for initial- 
condition errors applied simultaneously is equal to the sum of the results 
for each of these errors applied individually. 

These altitude errors are principally a function of time so that if 
a re-entry is of longer duration, the error will be worse. Figure 3("b) 
shows this effect for three different conditions. It should be noted 
that for both conditions of the standard vehicle the error at 300,000 
feet is around 1/3 of the total error at 100,000 feet, and that the error 
build-up occurs most rapidly during the last few hundred seconds. Thus 
the altitude error at the end of the run caused by incorrect initial 
conditions can be considerably reduced if the run can be made somewhat 
shorter, that is, if altitude can be determined from external means some- 
what sooner. 

The accelerometer biases and gyro drifts were studied separately 
and collectively, and in combination with initial- condition errors. In 
every case, the error due to a combination of error-producing sources 
was equal to the sum of the errors from the error sources individually . 



In addition, the error was directly proportional to the magnitude of the 
error source. Thus data which have been taken for one particular magni- 
tude for each error-producing source may be scaled according to the 
magnitude of the expected actual error source, and then the results for 
several error sources may be combined to give sui estimated value of total 
error . 

Figure 4 shows the altitude error versus time for different 
re-entries during which there was gyro drift. The reference time used 
in plotting this figure was that time at which the satellite passed 
through 300,000 feet altitude, which is where the air density effects 
first become noticeable. For those cases where = 2xl0~ 7 , there is 
no appreciable altitude error due to gyro drift until about 3OO seconds 
after t re f. This is to be expected, since gyro drift rotates the plat- 
form containing the accelerometers, and this has no affect until the 
accelerometers have an appreciable output (below 300,000 ft). The 
build-up of the altitude error due to gyro drift is a function of time 
of drift and the magnitude of the acceleration. For the 900,000-foot- 
altitude case, the re-entry time was considerably longer than for the 
677,000-foot case, which would imply greater errors. However, the 
acceleration level was, in general, considerably lower, implying smaller 
errors. These two effects combined to give a net reduction of error at 
any particular time following t re f but an increased error at the end 
of the run. 

The curve of figure 4 for the larger value of shows that the 

altitude error is proportional to the magnitude of the gyro drift . 

Figure 5 shows the altitude error caused by accelerometer biases. 

It can be seen that the altitude error caused by Uxp is smaller than 
that caused by U zp during the early portions of a run, but then builds 
up considerably larger in the later portions of the run. 

The effect of U zp on altitude error shows a reversal in sign. 

This is caused by the interaction of the bias, U zp , as it is resolved, 
the centrifugal acceleration term in the vertical equation (see eq. (11)), 
and the altitude error itself. For the longer runs, the bias error, as 
resolved through the range angle, and the resultant centrifugal accelera- 
tion error combine to cause the error reversal. For the shorter runs, 
these terms become insignificant much sooner, and the existing altitude 
error prevents the sign reversal. 

The same divergence is noticeable on the rims with Uxp 
accelerometer bias. However, it is not so pronounced since there is 
no sign reversal until the range angle becomes 180 . 

In an actual accelerometer, the quasi-steady component of error may 
vary slowly throughout a run, and will vary from rim to run. It is 
assumed, in figure 5, that this quasi-steady error is constant through- 
out a run, and is termed a bias. The value of bias chosen is the maximum 
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anticipated, steady-state error and, therefore, represents the worst 
possible case. Thus the curves for U^-p shown in figure 5, together 
with their images reflected about the zero error line (which would 
represent runs with negative bias) , represent an envelope of all possible 
errors for an accelerometer with an output uncertainty equal to the bias 
used . 

For errors in the platform Z accelerometer, U zp , the worst 
possible case is not represented by constant bias throughout the run 
because of the cancellation of error that occurs. However, until the 
range angle approaches jt/2 (approximately 1200 sec. for all vehicles 
and conditions) constant bias does represent the worst case, and so the 
curves on figure 5 for U zp bias, plus their images, do represent an 
envelope of possible errors up to 1200 seconds. Thus comparisons of 
errors caused by U zp should be made at a time prior to this range angle . 
The time at which the range angle becomes it/2 is indicated on this and 
subsequent figures. 

Figure 6 shows, for the standard vehicle and conditions, the effect 
on the altitude error of switching out the accelerometers for a while 
after t ; T4 is the time at which both the accelerometers were switched 
in again? Since t 3 = 27 seconds (for this value of AV) the curve with 
t 4 = 27 seconds is actually one for which the accelerometers are not 
switched out at all. Figure 6(a) shows that the error due to U^p bias 
can be reduced by switching both accelerometers out, and the optimum 
switch-in time appears to be at about ^00 seconds. From a comparison of 
this with the time histories in figures 2(a) and (c) and the altitude 
versus acceleration curves of figure 2(d) , it can be seen that at "J00 
seconds, f y is about half the bias level and f a is approximately equal 
to the bias level, so that if the accelerometers were turned on when their 
output reached the bias level, the results would be nearly optimum. 

Figure 6(b) shows the same thing for errors caused by U zp . Again 
the optimum switch-in time appears to be about ^00 seconds. 

Figure 7 shows the effect of turning off only one accelerometer 
when both and U zp errors are present, in comparison with curves 

where neither or both accelerometers are switched out. This comparison 
is valid only out to about 1200 seconds because of the change in sign of 
the error due to U zp . However, one can see that with both accelerometers 
switched out until t = 710 seconds the error at this point is smaller 
than one would get if either were switched out by itself. The case with 
both accelerometers in all the time appears to give smaller errors than 
with just f zm out. However, this is misleading, since, for the partic- 
ular values of bias chosen, there is a cancellation of the errors from 
U zp and Uxr). If one of the biases were negative, the errors would add, 
and the case with both accelerometers in would then give larger errors 
than the case with f zm out . 

Figure 8 shows the effect of switching out the accelerometers, and 
turning them on again individually when their resolved outputs (f a and f^) 



are twice the bias level. This is shown for a variety of conditions and, 
for each condition, is compared with the corresponding results when the 
accelerometers are not switched out at all. It can be seen that switch- 
ing the accelerometers out leads to reduced errors for all conditions . 

The time at which the range angle becomes n/2 is shown on each 
figure, and the errors from U Z p, beyond this point no longer represent 
the worst possible cases. 

Altitude errors at the end of the run (h = 100,000 ft), caused by 
various error producing sources acting singly, are summarized in table I 
for the standard vehicle starting from h Q = 677 ,000 feet with AY = 250 
feet per second. Also shown in the table are altitude rate error, range 
error, and range rate error. The errors caused by the accelerometer 
biases are shown with the accelerometers switched in all the time, and 
also switched out for the optimum time. Finally, the effect of all 
errors, simultaneously applied, for both conditions, is shown. Here, 
the errors are combined on an RMS basis, since one would not expect 
maximum errors from all sources simultaneously . 


CONCLUDING REMARKS 


From these data, it can be determined that an inertial navigation 
system can provide a re-entry vehicle with the required position and 
velocity information during the re-entry, in spite of the instability of 
the system equations . There will be errors but even for only fair iner- 
tial navigation equipment, with fairly good initial information, these 
errors do not become excessive. 

The errors from any source build up with time, so that those 
re-entries which take the least time will produce the smallest errors. 
The sooner radio contact is re-established with the ground the smaller 
the errors will be. 

Further, although switching out the accelerometers does reduce the 
error from accelerometer bias, the errors caused by other sources are as 
large as those caused by the accelerometers. Thus the complexity of 
switching out the accelerometers may not be warranted. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., March 21, I960 
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APPENDIX A 

DERIVATION OF ALTITUDE INSTABILITY 


The first of equations (5) and (10) can l* e used to show the 
altitude instability mentioned in the introduction. Since h « R r , 
the gravity term in each equation may be linearized, by means of the 
binomial expansion, as follows: 


Sr-Rr 


gr 


(R r + h)‘ 


1 + =- 


R- 


= Sr(l-f + 


(Al) 


Substituting this into equation (5) and a similar expression in hq 
into (10), and collecting terms in h or hj_ gives: 


Rr 


Sr 


h - h [ Q +2 ) = Rr9 - Sr “ ^7 


''i " (®i + ^ R r 


Rr®i " ®r - ^ 


zm 


(A2) 


When the vehicle speed is considerably less than orbital, so that 
9 2 « g r /R r , 1 these equations become: 


h - ^ h = -8 r - f 7 
R r 

hj_ - h^ = -g r - f zm 

i»r 


(A3) 


The equation for altitude error can be obtained by subtracting these two 
equations : 


2gr 

Rr- 


or 


•• 2g r 
h e - p~ 


h) - ”(fzm ~ 

f 7 ) 

(AU) 

h e = "(^zm “ 

f 7) 

(A5) 


1 At orbital speed = g r /R r 



When solved, this equation has an unstable root of the form 

t/ Rp/2g r _ j 

Ke with sJ~R T /2g T about 10 minutes. This period does not 

change appreciably until the vehicle velocity becomes about 0.8 of orbital 
velocity. 

This unstable solution for altitude error will have a magnitude K 
which is proportional to any initial altitude error, and to any subse- 
quent error in measuring acceleration. Thus one possibility for reducing 
the error caused by this instability is to reduce the accelerometer meas- 
urement error. In particular, if the anticipated measurement error is 
larger than the magnitude of the true acceleration, it would seem that 
zero would be the best value of acceleration to use in computing h^. 
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TABLE I.- ALTITUDE ERROR AT END OF ENTRY FOR A STANDARD VEHICLE 
ENTERING FROM 677,000 FEET WITH A AV OF 250 FEET PER SECOND 


Error at 100,000 feet 



Error source 



and 

magnitude 

ft 

9 e> 

radian 

V 

ft/sec 

0 e xio s , 
radian/ sec 

^eo 

= 500 ft 

5,427 

-0.0002 

10.1 

-0.094 

0 eo 

= 10” 4 radian 

-811 

.0001 

-2.9 

- .0292 

*60 

= 1 ft/sec 

1,598 

-.0001 

3-6 

-.00052 


mm Q 

= 5(10) radian/sec 

5,121 

-.0002 

8.9 

.1165 

“a 

_ y 

= 2(10) radian/sec 

-2,902 

-.0000 

-10.8 

.1207 

Accelerometer bias = 10 -4 g 





Uxp^ 

Accelerometers 

13,192 

-.0004 

23.9 

• 389 

u zp J 

t in all the time 1 

-883 

.0001 

•9 

.218 

^xp 

1 Accelerometers out 

6,583 

-.0002 

12.2 

.238 

U zp , 

[for optimum time 1 

112 

.0001 

1-3 

.098 

RMS sum 2 ''l 

Accelerometers in \ 

15,400 

.0005 

29.7 

.487 


all the time 
RMS sum 2 

Accelerometers out 
for optimum time 


10,400 


.0004 


1 The errors caused hy U Z p are given for a time of 1260 seconds 
(where Q = jt/2) . 

2 In computing the RMS sum, it was assumed that the errors caused 
■by u Z p remained constant subsequent to 1260 seconds . The actual 
error from the worst case was not determined and could he larger than 
the assumed error. 
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Figure 1.- Coordinates and vectors. 
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Figure 2.- Flight-path profile for different conditions and vehicles. 
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(b) Altitude versus time for other initial altitudes. 
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(d) Acceleration versus altitude. 
Figure 2.- Concluded. 
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Figure 4.- Effect of gyro drift on altitude errors for standard vehicle 
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Figure 5.- Altitude errors from accelerometer biases for different conditions and vehicles. 
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Figure 6.- Effect of accelerometer switch in time on altitude error for standard vehicle with basic 

conditions . 





Figure 6.- Concluded. 
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(f) Standard vehicle, h Q = 677,000 feet, AV = 475 feet per second 

Figure 8.- Continued. 
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Figure 8*- Concluded. 



